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Abstract 

The amplification of disk oscillations resulting from nonlinear resonant couplings 
between the oscillations and a disk deformation is examined. The disk is geometri- 
cally thin and general relativistic with a non-rotating central source. A Lagrangian 
formulation is adopted. The author examined the same problem a few years ago, but 
here we derive a general stability criterion in a more perspective way. Another distinct 
point from the previous work is that in addition to the case where the deformation 
is a warp, the case where the deformation is a one-armed pattern symmetric with 
respect to the equatorial plane is considered. The results obtained show that in ad- 
dition to the previous results that the inertial-acoustic mode and g-mode oscillations 
are amplified by horizontal resonance in warped disks, they also amplified by horizon- 
tal resonance in disks deformed by one-armed pattern symmetric with respect to the 
equatorial plane. If we consider local oscillations that are localized around boundaries 
of their propagation region, the resonance occurs at 4rg, where is the Schwarzschild 
radius. If nonlocal oscillations are considered, frequency ranges of oscillations where 
oscillations are resonantly amplified are specified. 

Key words: accretion, accretion disks — black holes — high-frequency quasi- 
periodic oscillations — relativity — stability — X-rays; stars 



1. Introduction 

Since the discovery of high-frequency quasi-periodic oscillations (QPOs) by RXTE satel- 
lite in neutron-star X-ray and black-hole X-ray binaries, much progress has been made both 
in observations of the high-frequency QPOs and in their theoretical modeling. This is because 
the QPOs are regarded to be a very powerful tool to explore the mass and spin of the central 
compact sources and also to explore the physical states of the innermost relativistic region of 
accretion disks. 

One of promising models of the high-frequency QPOs is that they are disk oscillations 
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resonantly excited in relativistic disks. At an early stage of observational development of the 
high-frequency QPOs, Abramowicz and Kluzniak (2001) and Kluzniak and Abramowicz (2001) 
already pointed out the importance of resonant processes as the cause of the high-frequency 
QPOs. Their model is that the high-frequency QPOs are the results of resonant couplings 
between the vertical and horizontal epicyclic oscillations in disks. This idea is very fascinating, 
but there seems to have still uncertainty concerning their excitation processes. 

Kato (2003, 2004), on the other hand, proposed another resonant model of the high- 
frequency QPOs, where the QPOs are regarded to be disk oscillations resonantly excited on 
deformed disks. In this model, a deformation of disks from an axially-symmetric equilibrium 
state is an essential ingredient for resonant excitation of disk oscillations. He especially consid- 
ered the case where the disk deformation is a warp. An outline of the model is as follows. A 
nonlinear coupling between a disk oscillation (hereafter we call it the original oscillation) and 
a warp brings about some disk oscillations (we call them intermediate oscillations). The in- 
termediate oscillations make a resonant coupling with the unperturbed disk at some particular 
radius of the disk. After this resonant coupling with the unperturbed disk, the intermediate 
oscillations feedback to the original oscillation by nonlinearly coupling with the warp. Since 
the nonlinear feedback processes involve a resonance, the original oscillation is amplified or 
dampened. For this nonlinear coupling to occur, however, some conditions must be satisfied. 
These conditions specify the place(s) where the resonance occurs and the frequencies of res- 
onant oscillations (Kato 2004). Among such resonant oscillations, some ones are amplified 
and others are dampened (Kato 2004). Application of these resonantly-excited oscillations to 
the high-frequency QPOs has been made by Kato and his collaborators, e.g., Kato and Fukue 
(2006) and Kato (2007). 

Examination in which cases the oscillations are resonantly excited was made by Kato 
(2003) and (2004), but their analyses are very lengthy and mathematically not perspective. 
The reason is that cylindrical coordinates are introduced at a early stage of formulation. In 
spite of this, the results obtained are rather simple. This suggests that a more perspective and 
general derivation of stability criterion is possible. Really, we find that the stability criterion 
can be derived in a simpler and perspective way. A purpose of this paper is to show this. The 
results obtained confirm the stability criteria obtained by Kato (2004). 

In addition, we show in this paper that a warp is not only one type of deformation that 
can resonantly amplify disk oscillations. One-armed deformations of disks that are symmetric 
with respect to the equatorial plane are also one of possible candidates that can resonantly 
amplify disk oscillations. 

In our previous work, we specifically consider the cases where oscillations are local and 
localized around the boundaries of their propagation region. In this paper we relax this ap- 
proximation, and specify the frequency ranges (in other words, the radial ranges) where global 
oscillations are resonantly amplified. 
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2. Summary of Basic Nonlinear Hydrodynamic Equations 



Here, we summarize basic nonlinear hydrodynamic equations that are necessary to ex- 
amine nonlinear resonant processes and their stability. Details are given by Kato (2004). 

We are interested in general relativistic disks, since in the resonant instability processes 
that we consider here the general relativity is essential. The essence of the effects of the general 
relativity, however, can be taken into account within the framework of a pseudo-Newtonian for- 
mulation using the gravitational potential introduced by Paczyhski and Wiita (1980), provided 
that the central object has no rotation. 

The formulation presented here are quite general until some approximations are intro- 
duced later. The unperturbed state of the systems considered here is a steady equilibrium state. 
In Lagrangian formulation, the hydrodynamic equation describing perturbations over the state 
can be written as (Lynden-Bell and Ostriker 1967) 



Dt^ ^(-V*--Vpj, (1) 

where ^ is a displacement vector associated with perturbations, and Dq/DI is the time derivative 
along an unperturbed flow, Uq, and is related to the Eulerian derivative, d/dt, by 

In equation (1), 6 represents the Lagrangian variation of the quantities in the subsequent 
parentheses, and ip is the gravitational potential. Other notations in equation (1) have their 
usual meanings. 

The basic equation (1) is valid even when the displacement ^ has a flnite amplitude. 
The purpose here is to explicitly write down the right-hand side of equation (1) up to the 
second-order quantities with respect to Now, we consider a system consisting of a central 
object and a disk surrounding the central object. Perturbations are only on the disk and the 
self-gravity of the disk gas is neglected. That is, the gravitational potential, ip{r), comes from 
the central object, and there is no Eulerian perturbation of iplr)} 

To represent the Lagrangian variation of density, 6p, by quantities up to the second- 
order ones with respect to ^, we need the equation of continuity. We assume here that the 
perturbations are adiabatic. Then, the Lagrangian variation of pressure, 6p, can be expressed 
by quantities up to the second-order ones with respect to Then, the nonlinear hydrodynamic 
equation describing adiabatic, Nonself-gravitating perturbations is written as, after lengthy 
manipulation, 

po-^ + 2po{uo • V)^ + m = PoC(|,|), (3) 



^ The Lagrangian formulation by Lyndcn-BcU and Ostriker (1967) is quite general, and the effects of self- 
gravity can be taken into account. In this paper, however, we neglect the effects of self-gravity, since we are 
interested in standard accretion disks. 



where L[^) is a linear operator with respect to $, and is (Lynden-Bell and Ostriker 1967) 

L{^) = Poiuo ■ V)(no • V)^ + po(l ■ V)(V^o) + V [(1 - T^)podiv^ 

-poV(div^) - V[(£ ■ V)po] + ■ V)(Vpo), (4) 

and po{r) and po{r) are the density and pressure in the unperturbed state, and Fi is the 
barotropic index specifying the hnear part of the relation between Lagrangian variations 6p 
and 6p, i.e., (5p/po)iincar = ri(<^p/Po)iincar- The right-hand side of wave equation (3) represents 
the nonlinear terms, which consist of two parts, i.e., p^C = poC^ + poCp. The former comes 
from —pq6{ViP) in equation (1) and poCp from —p^SiVp/ p). Here and hereafter, for simplicity, 
we consider the case of Fi = 1. Then, we have (Kato 2004) (an expression for nonlinear terms 
in the general case of Fi 7^ 1 is given in appendix) 
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PoCp = -^U|^Ve.). (6) 



and 

OTi \ or J 

3. Wave Energy and Work Done on Oscillations 

We now assume that the disks are deformed from an axisymmetric steady state by 
some external or internal cause. The Lagrangian displacement associated with the deforma- 
tion is denoted by ^^(r,t). The deformation may be time-periodic with frequency 00^ as 
= exp(iu;^t)^ (r). In order to avoid unnecessary complication, however, we take = 0, 
hereafter. The effects of non-zero can easily be taken into account in the final results. 
As the deformation, a warp will be most probable, but it is not only one candidate of possi- 
ble deformation. A plane-symmetric one-armed spiral deformation is one of another possible 
candidates that excite disk oscillations, as discussed later. 

Our purpose here is to examine how the behavior of disk oscillations are affected by 
the disk deformation. As shown below, some of oscillation modes are resonantly excited on 
deformed disks through nonlinear coupling with disk deformation. The nonlinear coupling 
processes are schematically shown in figure 1 of Kato (2004). 

Here we consider a disk oscillation mode. The displacement vector associated with the 
oscillation is denoted by $, with frequency u, i.e., ^ = exp{iujt)^{r) , in the limit of no deformation 

of disks. The first step of the nonlinear interaction between the disk oscillation characterized 

^ w 

by (cu,^) and the deformation characterized by (0,£ ) introduces two kinds of intermediate 
oscillations. Let us denote the displacement vector associated with these intermediate oscilla- 
tions by ^±'^[= exp{iut)$,_j_ ]. In other words, we characterize the two intermediate oscillations 
by {uj,^_^ ) and {uj,^_ ), where represents the intermediate oscillations resulting from the 
coupling between ^ and ^ , while ^_ does those resulting from the coupling between ^ and 



^ , where the asterisk represents the complex conjugate. That is, and ^_ are described, 
respectively, by 

-u'poC + 2tcopo{uo ■ V)C + lCC) = llPoCilt) + PoCU^M (7) 

-uj'pot + 2iu;po{uo ■ V)t + Lit) = l[PoC{lt*) + PoC{t*M (8) 

Next, the second stage of the nonlinear coupling is considered, which is a feedback process 
returning to the original oscillation, ^, by (or ^™*) interacting with The feedback is 
described by 

1 int W* W* int 

-cuVol + 2^a;poK-V)^ + X(|) = -[poC(£+ ,1 )+PoC(^ )], (9) 

- int " W 

in the case where the nonlinear coupling between and $, feedbacks to the original oscilla- 
tion, ^. On the other hand, in the case where ^^^^ and couples to feedback to the original 
oscillation, the equation corresponding to equation (9) is 

-c.Vo^ + 2^c.poK- V)^ + X(0 = l[poCCC\t) + PoC{t,C)]. (10) 

An important point to be noted here is that as a result of this feedback process, the 
original disk oscillation is amplified or dampened, i.e., the frequency, u, can be no longer real, 
since in the feedback a resonant process is involved as shown later. How much is the imaginary 
part of uj7 This can be examined from equations (9) and (10) by using the fact that the 
operators ipo('"o ■ V) and L are Hermitian (Lynden-Bell and Ostriker 1967). 

Let us now multifly equation (9) by ^ and integrate over the whole volume, and consider 
the imaginary part of the resulting equation. The volume integration of $, ■ L{$,) does not have 
any imaginary part, since the operator L is Hermitian. Now we write uj = ujq + iuji, where Ui is 
the imaginary part of u. Then, from equation (9) we have 

-2^i| Poe[u^o-^{uo■VMdV = Q^J Po^^*[CCC\t*) + cCt\C)W- (H) 
This equation has a clear physical meaning. The integral on the left-hand side is related to the 
wave energy. Really, the wave energy is expressed as (Kato 2001) 

E = ^uoJ Po^^*[uo-i{u-V)]idV. (12) 

The right-hand side of equation (11) is thus —4{uJi/u!Q)E. Hereafter, the subscript to u; is 
neglected for simplicity. In the case of oscillations that have m-arms in the azimuthal direction 
in geometrically thin disks, the wave energy E can be expressed in terms of $ as (Kato 2001) 

E = ^Juj{uj- mn)pom^ + \i\'')dV, (13) 

where the cylindrical coordinates (r, if, z) whose origin is at the disk center and the 2-axis is 
the axis of disk rotation have been introduced. 
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Next, the right-hand side of equation (11) is considered. This is related to the work 
done on oscillations. Let us assume that an external force / on the oscillations acts per unit 
volume. Then, the work done by the force on the oscillations per unit time, W, is 

W={Jpo^-fdV), (14) 

where ( ) represents the time average. If we write ^ = 3fJ[exp(zu;t)|] and / = 3?[exp(zu;t)/], we 
find that W can be expressed as 

^ = Y^j Pol* ■ fdV. (15) 

This means that in the present case, i.e., equation (9), the work done on the original oscillations 
by the nonlinear coupling terms is 

UJ C \ * ^ int ^ W* ^ W* ^ int 

W^ = --sj'-p,i[C{i^,i ) + C{i ,i^)]dV, (16) 

where the subscript + is attached to W in order to emphasize that the feedback occurs through 
. Hence, equation (11) can be expressed as 

= (17) 

which is an expected result. 

^ int 

In the case of the coupling through ^_ , the work done the original oscillations is written 

as 

w- = iporic(t°*.r) + ci£',t)\iv, (18) 

and 

w 

= (19) 

The quantitities cui's in equations (17) and (19) are different by difference of coupling path. In 
order to distinguish them, the subscript + or — is attached to ui. 

The purpose in the following sections is thus to express W± more explicitly and to 
evaluate the sigh of W±. It is noted that there are two cases where the wave energy is positive 
and negative. If the waves exist dominantly inside the corotation radius (i.e., the radius of 
UJ = mQ), the wave energy E is negative [see equation (13)] and a negative work (W < 0) leads 
to instability {—Ui > 0), while the waves exist mainly outside the corotation radius, we have 
E > 0, and W > leads to instabihty. 

The coupling term C consists of two terms of and Cp, given by equation (5) and 
(6), respectively. From equaion (5) we easily find that 

/ Poi*c,cc\t*)dv = I poCc,a*r^"*)dv. (20) 

In the case of the integration of po^ Cp over the volume, we perform two times the integration 
by part. Then, we have 
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/ P,i'cAS':\t")dv = / p^tcAi'.t'W. (21) 

Hence, W+ given by equation (16) is reduced to 

w.=^^j ip„«:V(r.r*) + ccr.tm (22) 



^ * int ^ W 

It is very important to note here that ^ , ^ , and ^ in integrations (20) and (21) are comm- 
nesurable each other as one of their examples is shown above. This nature is not hmited in 
the case of Fi = 1. This holds generally even in the case of Fi 7^ 1, as shown in appendix. The 
terms in the brackets of equation (22) are just the complex conjugate of the right-hand side of 
equation (7), which makes the calculation of W+ easy, as shown later. 

In the case of the feedback through ^_ , we similarly have, from equation (18), 

UJ r 1 ^int ^*^W ^W"* 

W. = -qJ-Po^^[C{^,^ ) + C(| ,t)]dV. (23) 

The terms in the brackets are just the complex conjugate of the right-hand side of equation 
(8), which makes the calculation of W- easy, as in the case of W+. 

4. Dispersion Relation 

Hereafter we restrict our attention to oscillations in geometrically thin disks. The steady 
unperturbed disks are assumed to be axially-symmetric and to have no motion except for rota- 
tion. Here, cylindrical coordinates (r, ip, z) are employed, in which the z-axis is perpendicular 
to the disk plane and the origin of the coordinates is at the disk center. The rotation is then 
described as Uq = {0,rQ{r),0). When we need numerical figures, we adopt the Keplerian form 
of Q, i.e., Q = Qk, since geometrically thin disks are considered. 

The oscillations that we consider here are assumed to have m arms in the azimuthal 
direction, i.e., ^ oc exp(— zmy?). We further assume that the radial wavelength of the oscillations 
is moderately short, and the local approximation in the radial direction can be adopted. That 
is, ^ oc exp{ikr) and the radial variations of physical quantities in the unperturbed state are 
neglected compared with the radial variation of ^. Concerning the vertical direction, however, 
we cannot adopt local approximations, since the disk is geometrically thin. 

For mathematical simplicity, we assume that the disks are vertically isothermal. In the 
vertically isothermal disks the density pQ{r,z) is stratified as (e.g., Kato et al. 1998) 



po{r,z) = poo[r)exp 



2iS;)]- 

where poo is the density on the equatorial plane, and H is the half-thickness of the disk and is 
related to the vertical epi cyclic frequency, Q±, by 

nlH^ = ^ = cl{r). (25) 
Po 

The vertical epicyclic frequency, Q±, is equal to the angular velocity of the Keplerian rotation, 
f^K, in the case of the central object being non-rotating, and practically equal to Q, since the 
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disk is assumed to be geometrically thin. Hereafter, however, we use without using or 
Q, so that we can trace back the effects of Q± on the final results. Furthermore, we assume 
that oscillations also occur isothermally. Then, we can neglect the term of V[(l — ri)podiv^] 
in the operator L in equation (4). Under these approximations, we express the r-, ip-, and z- 
components of the homogeneous parts of wave equation (3). Then, we have 

[-{uj - mfi)2 + K^- + k'^cl]ir - i2n{uj - mn% - ik(^cl^ - zn\i)j = 0, (26) 

-{uj-mn)'^i^ + i2n{uj-mn)ir = 0, (27) 

-(. - mn)% - ^kcl^- [cf^ - nlz^j^ - file;) = 0. (28) 

These wave equations are differential equations, and it is troublesome to rigorously 
solve these equations. In vertically isothermal disks, however, the wave equations can be 
approximately solved by separating r and z dependences of ^'s as, for example, ^^ir^z) = 
f{r)g{z/H), where / and g are some functions of r and z/H, respectively (Okazaki et al. 
1987). Here, the functional form of g{z/H) is determined by solving an eigenvalue problem 
describing the vertical behavior of the oscillations, and we have 



g(z)^nn^-j, (29) 

where ?i„ is the Hermite polynomial of argument z/H, and n{= 0,1,2,...) characterizes the 
number of node(s) of oscillations in the vertical direction. More explicitly, we have 

Ur)=l^^{r,^)nn{z/H) (30) 

Ur)=Unir,v)nn{z/H), (31) 

Ur)=iAr,v)'>^n-i{z/H). (32) 



It is noted that the number of node(s) of C,z in the vertical direction is smaller than those of 
and by one, as shown in equation (32). However, the subscript n (not n — 1) is attached 
to as ^z,n in order to emphasize that ^r,n, ^cp,n, and ^z,n are a set of solutions. The fact that 
the ^-dependences of ^r, ^if, and given by equations (30) - (32) are really solutions of the 
set of equations (26) - (28) is found by substituting equations (30) - (32) into equations (26) - 
(28) and conforming that the right-hand side of equations (26) - (28) can be packed into terms 



proportional to Tin, and 7i„_i, respectively. The results are (Kato 2004) 

[-{u - mQf + 4fi2 + A;2c2]4n - ^21](cj - ml])4,n + i{kH)n\l^^n = 0, (33) 

— (u — miiy^^n + i2^l{uj — mf2)^r,n = 0, (34) 

[-{iu - mnf + nnl%,n - in{kH)nl^r,n = 0. (35) 
The solvability condition of these homogeneous equations is 

D{uj,m,n) = [{lu - mnf - k^] [{uj - mVlf - nn\] - k'^cl{uj - mVlf = 0. (36) 



This is the well-known local dispersion relation of oscillations characterized by uj, m, and n. 
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5. Expressions for Work Done on Oscillations and Stability Conditions 



Our purpose here is to evaluate W± given by equations (22) and (23). As disk deforma- 
tion, we consider two cases: 

• case (i): rn^ = 1 and ri^ = 1. 

• case (ii): rri^ = 1 and = 0. 

The former is a warp and the latter is a one-armed pattern plane-symmetric with respect to 
the equatorial plane. 

In both cases, for simplicity, the radial wavelength of the deformation is assumed to 

be rather long compared with the radial wavelength of the disk oscillations. That is, the 

^ w 

radial variation of ^ is neglected compared with that of ^. For convenience, the part of the 
2;-dependence of ^ is separated from other parts as 

eW(r)=e;(r,vp)Kw(^/i/), (37) 

i:^{r)=Ur.mn-^iz/H), (38) 

iTir) = Ur.v)^n-'-i{z/H), (39) 

where ri^ = 1 in case (i), and ti^ = and is case (ii). 

5.1. Inhomogeneous Wave Equation for Intermediate Oscillations 

Nonlinear coupling between ^ and introduces intermediate oscillations, represented 
by = exp(za;t)^ . Related to the fact that TiiKn = T^n+i + n'Hn-i, there are two major 
normal modes of intermediate oscillations induced by the coupling between oscillations and a 
warp. They are those characterized by Tin+iiz/H) and those characterized by T-Cn-i{z / H) ."^ 
Concerning the (/^-dependence of $, , there are also two cases where $, is proportional to 
exp[— z(m + !)(/>] and to exp[— (m — l)(p]. That is, there are four modes of ^ by differences of 
the ip- and ^-dependences in the case of a warp. In order to treat these four modes separately, 

'^int 

4 is written as 

t\r) = Y.C%dr,m^iz/H), (40) 

ii\r) = T.^l^X,dr,mf^i^/H), (41) 

t\r) = ECl,.(^'^)^^-i(^/^)' (42) 

where the subscript ± is attached to to represent the two cases where ^ is proportional 
to exp[— i(m ± 1)(/)] (see also section 3), and the subscript n represents n + 1 or n — 1 in case (i) 
and n = n in case (ii). 

Now we solve inhomogeneous wave equation (7) of . All the terms on the right-hand 
side of equation (7) have the (/^-dependence of exp[—i{m + l)(p]. Their ^r-dependences, however, 



^ In case (ii), the z-dependence of the intermediate oscillations is Tin{z/H) alone. 
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are not always limited to T-Cn{z/H) nor T-ln^i{z/H), although the major terms are those with 
such ^-dependences. Here, we write the right-hand side of equation (7) as 

+ C(r,0]r = PoE4+,"(^'^)^n(^/^) + ■ ■ • (43) 
^ n 

IpoiCCtt) + C(r,0]^ = Poj:^^,+Ar,V^)nn{z/H) + ■■■ (44) 

n 

IpolCCtt) + C{tX)]z = p,Y.K+Ar.mn-i{z/H) + ■ ■ (45) 

n 

and concentrate our attention on the terms shown by A's. Here, + ■ ■■ denotes other terms 
orthogonal both to Tin and Ha-i, and the subscript -|- is added to A's in order to emphasize 
that they have the (/^-dependence of exp[— 'i(m + l)(p]. 

Under these preparations, we obtain from equation (7) equations describing ir,+,h, ^ip,+,h, 
and ^z,+,h as [cf., equations (33) - (35)] 

{-[(.;- (m + + k'- An'+k'clY^ln - - ("^ + l)^]CV,n 



+tikH)n'J::l^^ = A,+,^, (46) 

-[lo- {m + + z2n[u; - (m + l^^^l^^ = (47) 

{-[(.;- (m + + nl^ljei^,, - zn{kH)Ql^^l^^ = (48) 
where n = n + lorn — lin case (i), and n = n in case (ii). 



^ int 

In the case of intermediate oscillations characterized by ^_ , we can derive similar equa- 
tions as the above. That is, we denote the right-hand side of equation (8) by A„ with similar 
expressions as equations (43) - (45). Then, as equations describing ^r-,h, ^(/3,-,n,and ^z-,n, we 
obtain the same equations as (46) - (48), respectively, except that m -|- 1 is changed to m — 1 
and the subscript -|- in all variables are now changed to the subscript — . 

The equations describing the intermediate oscillations, i.e., equations (46) - (48), and 
the corresponding equations for ^_ - can easily be solved with respect to ^T-,±,ri; ^ip,±,h, and 
ir,±,n- The results show that they are expressed in a form of N±/D±, where A^+, for example, 
is some combination of A_|_ fj's on the right-hand side of equations (46) - (48) and D± is 

D±(u;,m±l,n) =|-[cu - (m ± 1)^]^ + K^jj - [cu - (m ± + nfi^j 

~[uj-{m±l)n]''k^cl (49) 

5.2. Resonance 

The appearance of D± in the denominator comes from the fact that the homogeneous 
wave equation for the intermediate oscillations gives the dispersion relation as the solvability 
condition and this is D±{uj,m±l,n) = 0. In other word, at the radius where D±{uj,m±l,n) = 
is satisfied, the intermediate oscillations resonantly response to the external forces given by the 
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right-hand side of the wave equations. In the case of intermediate oscillations of the 
external forcing terms are the right-hand sides of equations (46) - (48). 

Since we are considering the cases where k'^c^ is much smaller than Q'^ (i.e., geometrically 
thin disks), the last term of the expression for D± can be approximately neglected (particle 
approximation). This means that the resonance occurs at two radii. One occurs at a radius 
close to the radius where 

[cu — (m ± l)f2]^ — = (horizontal resonance) (50) 

holds, and the other does at a radius where 

[u — {m ± 1)Q]'^ = hfl^. (vertical resonance) (51) 

The former is called hereafter the horizontal resonance since the resonance occurs by horizontal 
motions, and the latter the vertical resonance since it occurs by vertical motions. Our purpose 
here is to evaluate the right-hand side of equation (22) and (23) in these two cases. In order to 
explicitly see the effects of resonance, we introduce variables Cr,±,n, C<f,±,h, and Cz,±,h defined by 

Then, considering that (l/2)po[C'(^ ,^ ) + C(^ ,^ )] is just the complex conjugate of the 
right-hand side of inhomogeneous wave equation with respect to [see equation (7)], and 
that the inhomogeneous wave equations are explicitly written by equations (46) - (48), we have 

;M+-[c(i,r )+c(r ,1)] 

- [CU - (m + l)fi]2 + _ 4fi2 + eclYr,+,nQ- 

-Z2Q[U; - (m + l)n]{(r,+,nC;,+,n - C,+,n4,- 
-ink Cg ( — ^r,+,hCz,+,h ~^ Cr,+,nCz,+,h) 



2' 



(53) 



Our purpose here is to evaluate the imaginary part of the volume integration of the right- 
hand side of equation (53) to have an expression for W+. As we see later, D± is proportional 
to r — Tc near the resonant radius of D± = 0. Here, is a function of u, m, n. The terms 
in the brackets of equation (53) are all real at a glance. Hence, the radial integration of the 
right-hand side of equation (53) has no imaginary part, unless the terms in the brackets have 
terms proportional to D^. If they have terms proportional to D^, however, the right-hand 
side of equation (53) has terms proportional to 1/-D+ with real coefficients. Since, 1/-D+ is a 
pole in the complex r-plane, the radial integration of 1/-D+ brings about an imaginary part, as 
carefully discussed later. 

In the case of the horizontal resonance, we have 
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-[uj-{m+ mf + /.^ + k\l ~ (54) 
while in the case of the vertical resonance, we obtain 

Under these considerations, we pick up only, from the right-hand side of equation (53), the 
terms which contribute to the imaginary part when they are integrated in the radial direction. 
Then, we have 

Wn,+= --i j -p,i^ [Cii ,i ) + c{i ,i)]dv 

^/o ^3/2~,c- [ ^ 1 2 rpoo{r)H 

for the horizontal resonance. The final equality is obtained by performing the vertical integra- 
tion by using the nature of the Hermite polynomials. Here, the subscript H has been attached 
to and in order to emphasize that the horizontal resonance is now considered. In the 
case of the vertical resonance, we have 

LJ f ^int ^*^W ^W^* 

w^v,+= |s / M+ [c(^ ,i^) + c{i^,i )]dv 

.|„.)3.(«-i):./^lcV.r^^*, (57) 

where the subscript V has been attached to Wj^ and in order to clearly show that the 
vertical resonance is considered. It is noted that in the case of the vertical resonance, W , is 
proportional to (n — 1)!, not h\. 

In the case of coupling through ^_ , we have similar expressions as the above. They are 

W^H,-^|(2.)3/^.i./^|6..P3^^^ (58) 

W^v,- = ^{2^r\n - m I -T^\C.-,^-p^dr. (59) 



2' ' ' 'J Dv,- fi:2-n(]2 
It is noted that n = n + l orn — 1 in the case where the disk deformation is a warp [case 
(i)], while fi = n in the spiral deformation of disks which are symmetric with respect to the 
equatorial plane [case (ii)]. 

5. 3. Integrations of 1/D around the Resonant Radius 

The imaginary part of equations (56) - (59) comes from the integration around the 
resonant radius, since the radius is a pole of the integrant in the complex r-plane. To perform 
the integration around a pole, we must be careful. If the frequency of oscillations, uj, is taken 
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to be real, the pole is just on the real axis of r, and it is a problem what path we should 
take around the pole, when we make integration along the real r-axis. It is known in plasma 
physics, however, that the following procedure can satisfy the requirement of causality. That 
is, we perform the integration along the real r-axis, tentatively assuming that the oscillations 
in consideration are growing (i.e., cj; < in the present problem), and then the results obtained 
in this way are analytically extended to the case of cj; > 0. Hereafter, we examine the horizontal 
and vertical resonance separately with this standard procedure, 
(a) Horizontal resonance 

For the horizontal resonance, we have [u; — (m ± — ~ and around the resonant 
radius, Vc, the first term of the Tayor expansion of -Dh,± gives 



Dn,± = 2{mi-K')Gu,± 



r — rr 



.UJ - 



(m±l)a 



where 



Gh,± =|(m ± l)n[u - (m ± i)n] 



dlnQ 
dlnr 



2 dlnn I 
dlnr J c 



(60) 



(61) 



and the subscript c denotes the values at the resonant radius. Considering these expressions 
and assuming <0 (the requirement of causality), we have (Kato 2004) 
dr TTTr 



D 



H,± 



'2(nfi2 



iSign[ti; — (m ± l)^l]c- 



(62) 



It is noted that this result is valid even when cui > as mentioned before, although it is obtained 
under the assumption of c^i < 0. 

Combining equations (56) [or (58)] and (62), we have, for horizontal resonance. 



--(27r)3/2vrn! 
4 



-sign[co' — (m ± l)^]c- 



(63) 



This result shows that the sign of the work done on the oscillations by the horizontal resonance 
is determined by the sign of u — {m ± 1)Q at the resonant radius. The work is positive if 
UJ — {m ± 1)Q at the resonant radius is negative, 
(b) Vertical resonance 

The work done on the oscillations by the vertical resonance is obtained by a similar way 
discussed above. That is, in the case of vertical resonance, the resonance occurs at a radius 
close to [cj — (m ± 1)^2]^ — fiQ^ = 0. Around this resonant radius, rc, we have 

'r — rr 



Dv,± = -2(nfil-/€')Gv, 



^ uji 



G 



v,± 



where 



Gy ± =|(m ± - (m ± 1)1^]^ + 2n( ] . 

I dlnr V amr / J c 



(64) 



(65) 



The radial integration of 1/ -Dv.± along the real r-axis then gives (Kato 2004) 
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dr 



rsign[u; — (m ± l)Q]c 



Then, from equations (57) [or (59)] and (66), we have, for the vertical resonance. 



--(27r)3/2 



4 



IV. 



(nfi2 



IC 



2,±,n I 



v,±| 



-sign[ci; — (m ± l)f2]c 



(66) 



(67) 



It is noted that in the vertical resonance the sign of work done on oscillations is also determined 
by the sigh of — (m± 1)1] at the resonant radius. The radius of resonance is, however, different 
from that in the case of the horizontal resonance, although both radii of the horizontal and 
vertical resonances are simply denoted by the same symbol Tc, without distinguishing them. 

5.4- Stability Condition 

The wave energy given by equation (13) can be expressed as 
(27r)3/2 



-uj{uj -m^l)c{r Hpoo)cEr, 



where En is a dimensionless quantity given by 
rHpQQ uj — niQ 



E„ 



n! — — + n 



(68) 



(69) 



(riJpoo)c - mil)c \ 

Hence, the growth rate, — tui [see equation (17) or (19)], can be expressed in the case of the 
horizontal resonance as 

{Cr,±,n/rf 1 sign[u;-(m±l)fi]c 



-'^i.H.zt = --n\ 



E„ 



{uj-mQ){nnl-n^y\Gn,±\ 
Similarly, in the case of vertical resonance, we have from equation (67) 

iCz,±,n/rf ] sign[t^- (m± l)fi]c 



4^ 



E„ 



(70) 



(71) 



{uj-mn){nni-K'^y\Gn,±\ 
The growth rates given by equations (70) and (71) show that the resonance leads to growth of 
oscillations when the sign of [u — {m±l)Q]/ {u — mQ) at the resonant radius is negative. In the 
foUowings, we examine in what types of resonance and in what wave modes the sign becomes 
negative. 

6. Resonant Radii and Resonant Frequencies of Locally Unstable Oscillations 

This examination is separately made below for various combinations of types of resonance 
and of wave mode. 

6.1. Horizontal Resonance of Inertial- Acoustic or G-Mode Oscillations 

In the horizontal resonance, the resonance is realized at the radius where [u; — (m±l)f2]2 — 
~ [see equation (50)]. The inertial-acoustic or g-mode oscillations will be predominantly 
exist at the radius where {uo — mQ)"^ — ~ is satisfied (see e.g., Kato and Fukue 2006). For 
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resonance to occur efficiently, both radii of [uj — [m ± = and {uj — mfi)^ = must be 
the same. This condition is satisfied at the radius of 

^ = ^ (72) 
in two cases where 

case (a) : u — {m + l)Q = —K, and uJ — mQ = K 

case (h) : uj — {m — 1)Q = K, and 00 — 771^ = —^. (73) 

We see that in both cases the sign of [u — {m + {u — niQ) or [u — {m — {u — mQ) is 
negative. This imphes that —00^ > and the resonance amphfies the oscillations. 

Let us see more in detail. In the former case of (a), [uj — {m + l)fi]c < and the work 
done on the oscillations is positive [see equation (63)], and the wave energy is positive [see 
equation (68)]. Hence, the oscillations grow. In the latter case of (b), [u — {m — l)Q]c > and a 
negative work is done on the oscillations. The oscillations, however, grow since the wave energy 
is negative. Various quantities related to resonance and amplification are summarized in tables 
1 and 2. 

The above argument can be applied not only to the case where the disk is warped (i.e., 
171^ = 1 and = 1), but also to the case where the disk is deformed by a one-armed pattern 
which is symmetric with respect to the equatorial plane (i.e., rri^ = 1 and ri^ = 0). 

6.2. Horizontal Resonance of C-Mode or Vertical P-Mode Oscillations 

These oscillation modes exist predominantly at the radius where [u — mfi)^ — nVL]_ ~ 0. 
Both conditions of [u; — (m± l)f2]^ = and (a; — mfi)^ = nVt\ can be simultaneously satisfied 
at the radius of 

K = n^/^n± - n. (74) 

when 

case (a) : u; — (m + l)f2 = K and u — niQ = n^^'^Q± 

case {h) : UJ — {m — 1)Q = — hi and uj — niQ = —n^^'^Q±. (75) 

This consideration shows that at the resonant point the sign of [u — {m ±l)Q]/{uj — mQ) is 
positive, showing that — cui^n.i given by equation (70) is negative. That is, the resonance 
dampens the oscillations (see tables 1 and 2). 

6. 3. Vertical Resonance of Inertial- Acoustic or G-Mode Oscillations 

The condition of the vertical resonance is [cu — (m ± l)f2]^ ~ hVt\. The radius where the 
inertial-acoustic or g-mode oscillations predominantly exists is {uj — mQ)"^ = k^. Combination 
of these two relations shows that the resonance occurs at 

K = h^/^n^ - n, (76) 
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when 

case (a) : — (m — 1)^2 = n^^'^Q± and u — mQ = k 

case (b) : — (m + = — n"^/^f2_L and uj — 'mQ = —K. (77) 

In these two cases the sigh of [u — (m ±1)1]]/ {u — mQ) is positive. Hence, equations (71) shows 
that the resonance dampens the oscillations. 

The vertical resonance of c- and vertical p-mode oscillations is absent, since the condi- 
tions of (cu — m^iy — h^^'^^l± = and (cu — mf2)^ — n^/^fij, = can not be simultaneously satisfied 
at any radius. 

6.4- Rough Estimate of Growth Rate 

Finally, the order of growth rate is briefly estimated in the case of horizontal resonance 
of inertial-acoustic or g-mode oscillations in warped disks. To do so, we must evaluate the 
order of The nonlinear coupling terms between ^ and are given by equations (5) 

and (6). Some detailed calculations show that the radial component of poCp, for example, 
is on the order of pQ^^-riVL^ [kH) / (Kato 2004), where and r] are the orders of the radial 
displacement associated with the oscillations and of the vertical displacement associated with 
the warp, respectively, and 1/k and are the radial wavelengths of the disk oscillation and 
the warp, respectively. Then, the wave equations describing the intermediate oscillations give 
the order of S,^^^ as 

0{t') = (78) 
Then, we have 

0{C')=CrV^'^. (79) 
The dimensionless wave energy contained in the disk. En, is found to be [see equation (69)] 

0(E.,-(^'^, (80) 

where L is the radial extend where the oscillations exist. Using these relations, we obtain from 
equation (70) that the growth rate is on the order of 

0{-Ui) = aikHyn, (81) 

where a = (?7/A^)^(rc/L), confirming the results obtained by Kato (2004).^ 

7. Resonant Excitation of Global Oscillations 

So far, we have considered resonant excitation of local oscillations, assuming that the 
oscillations are local in the radial direction in the sense that their radial wavelength is shorter 



In Kato (2004), is taken to be 
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Table 1. Resonant oscillations in warped disks (rn^ = 1, = 1). 



Incrtial-acoustic mode and g-mode oscillations 



type of resonance 


type of coupling resonant radius work 


wave energy 


stability 


horizontal 


m m + 1 — > m ^ /^e. = 4.0 positive 
m ^ 771 — 1 — > TO T /fg = 4.0 negative 


positive 
negative 


growth 
growth 


vertical 


TO ^ m + 1 ^ TO ?'/''g = 3.62, 6.46 negative 
TO ^ TO — 1 ^ TO r/rg = 3.62, 6.46 positive 


positive 
negative 


damping 
damping 


c-mode and vertical p-niode oscillations 


type of resonance 


type of coupling resonant radius work 


wave energy 


stability 


horizontal 


TO ^ 771 + 1 — > TO r/rg = 3.62, 6.46 negative 
TO ^ TO — 1 ^ TO r/rg = 3.62, 6.46 positive 


positive 
negative 


damping 
damping 


vertical 


absence 






Table 2. Resonant oscillations deformed by one-armed perturbations (i 


= 1, „w = 0). 




Incrtial-acoustic mode and g-mode oscillations 


type of resonance 


type of coupling resonant radius work 


wave energy 


stability 


horizontal 


TO ^ TO + 1 ^ TO r/rg_= 4.0 positive 
TO— i'TO — 1— >m fl^g= 4.0 negative 


positive 
negative 


growth 
growth 


vertical 


absence 






c-mode and vertical p-mode oscillations 


type of resonance 


type of coupling resonant radius work 


wave energy 


stability 


horizontal 


absence 






vertical 


absence 







than the characteristic radial scale-length of the disk (i.e., kr > 1). Certainly, we have exten- 
sively used the approximation in section 5 when the work done on oscillations is calculated. 
The processes of the analyses, however, suggest that the final result (i.e., the sign of the work 
done on oscillations by resonance is governed by the sign of u; — (m± at the resonant radius 
[see equations (63) and (67)]) is free from the approximation, although the magnitude of the 
work depend on the radial form of oscillations. This anticipation seems to be supported by 
considering the other limit of the particle approximation {kr > 1), although it is not proved. 

Different from the above, the resonant radii and resonant frequencies estimated in section 
6 come from the approximation that the oscillations are local. ^ Thus, in this section we relax 
the approximation and examine how the results in section 6 are modified in cases of global 
oscillations. 

In section 6 we have assumed that the oscillations exist predominantly around the boundary of the propa- 
gation region, since the group velocity of the local oscillations vanishes there. 
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We restrict, however, our attention only to the case where oscillations are inertial- 
acoustic or g-mode oscillations and resonance is horizontal, since we can show that in other 
cases oscillations are damped or the resonance appears only in the evanescent region of the 
oscillations. 

7.1. Horizontal Resonance of Inertial- Acoustic Oscillations 

Let us consider inertial-acoustic oscillations of frequency uj with azimuthal wavenumber 
m. The propagation region of such inertial-acoustic oscillations is specified by [uj — mQy > k^. 
In other words, the evanescent region of the oscillations is [u — mfi)^ < k^, i.e., mVL — k, < uj < 
mQ + K. The evanescent regions in the u-r plane (i.e., the so-called propagation diagram) are 
shown by dark area in figures 1 and 2 for m = 1 for m = 2, respectively. The region in the case 
of m = is supplement arily shown in figure 3. 

Next, we consider where the horizontal resonance occurs on the u-r plane. In the 
coupling through m — m + 1 ^ m, the resonance occurs on the curve of [u — {m+ 1)^^]^ = k^, 
while it occurs on the curve of [u — {m — = in the coupling through m ^ m — 1 ^ m 
[see subsection 5.2]. These curves are shown in figures 1 to 3 by thick (solid or dashed) or less 
thick dashed curves (see below for distinction of these curves). 

The first problem to be examined is in which part of the curve of [u — {m±l)Q]'^ = k^, 
resonance excites or dampens the oscillations. This can be examined by studying the sign of 
[u — {m± l)Q]/{u! —mQ) on the curve [see that W^s are given by equations (63) and (67) and 
the wave energy is given by equation (13)]. Above the curve of = mQ on the u-r plane, the 
wave energy, E, is positive, while it is negative below the curve. Next, in the case of coupling 
through m — > m + 1 — * m, the sign of the work done on oscillations is governed by the sign 
of a; — (m + In the case of coupling through m — m — 1 — » m, on the other hand, the 
sign of the work is governed by the sigh of u; — (m — Considering these situations, we find 
that resonance excites (amplifies) oscillations in the case where it occurs on the thick parts 
(including thick dashed parts) of the curves of [cu — (m± 1)^1]^ = k,^. However, resonance leads 
to damping of oscillations in the case where the resonant point is on less thick dashed parts of 
the curves. These results are shown in figures 1 and 2 for m = 1 and m = 2, respectively, and 
supplementarily in figure 3 for m = 0. 

Even when resonance excites oscillations, a strong excitation of the oscillations will 
not be expected in the case where the resonant point occurs in the evanescent region of the 
oscillations. This is because the coupling between the oscillations and the disk deformation is 
weak due to small amplitude of oscillations in the evanescent region. Hence, the parts of the 
curves of [u; — (m± 1)^2]^ = where resonance excites oscillations but the resonance occurs in 
the evanescent region of oscillations are shown by thick and dashed curves. 

As shown in figures 1 to 3, the thick solid parts of the resonant curves are only in a finite 
range of frequency. Propagations of some typical growing oscillations are shown in figures by 
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horizontal dashed hnes. The resonant radius of these oscillations is shown by a large open circle 
on the lines. It is noted that the cases considered in section 6 under the local approximations 
are those where resonance is assumed to occur at the crossing point of the resonant curve and 
the boundary curve of the propagation region, i.e., at 4rg. 

The oscillations labeled by cjll in figure 1 and by cul in figure 2 are trapped in a finite 
range in radius, bounded between the inner edge of the disk (~ 3rg) and the outer barrier of 
uj = Q — hi (for cjll) or of C(J = 2Q — n (for cjl). Compared with these oscillations, the oscillations 
labeled by (figure 1) and ci;hh (figure 2) are not bounded in their outer boundary. In this 
sense, the latter two oscillations will be less observable with definite frequencies, compared with 
the former two ones. This consideration suggests that the pair QPOs observed in some X-ray 
binaries are the oscillations labeled by cull and cul (see, e.g., Kato and Fukue 2006). 

It is noted that in the limit of local perturbations, the oscillations labeled by cull (figure 
1) and (figure 2) tend to the oscillations of case (b) of equation (73). Similarly, the oscillations 
labeled by c^h (figure 1) and whh (figure 2) tend to the oscillations of case (a) of equation (73). 

Finally, resonant excitation of axisymmmetric (m = 0) inertial-acoustic oscillations is 
briefiy mentioned.^ The results are shown in figure 3. In the case of m = 0, the resonance curve 
on the u-r plane is u = Q ± k, and the evanescent region of the oscillations is the dark area 
below the curve of K{r). The same considerations as those in cases of m = 1 and m = 2 show 
that axisymmetric inertial-acoustic oscillations are excited in the case where resonance occurs 
on the curve of u = Q — n. (On the thick dashed part of the curve, resonance occurs only in 
the evanescent region of the oscillations). The resonance on the curve of u = Q + k dampens 
oscillations (less thick dashed curve). The oscillations that can be excited by resonance are not 
trapped in a finite region. On the other hand, the oscillations whose frequencies are less than 
f^max can be trapped, bounded inside by the inner edge of the disk and outside by the barrier 
of uj = K. Such oscillations, however, do not have resonance in their propagation region. 

7.2. Horizontal Resonance of G-mode Oscillations 

The frequency-radius relation where resonance occurs for g-mode oscillations with fre- 
quency u and azimuthal wavenumber m is the same as that for inertial-acoustic oscillations 
with the same u and m. On the curve of the frequency-radius relation, the part where resonance 
dampens g-mode oscillations is also the same as that in the case of inertial-acoustic oscillations. 

A difference between inertial-acoustic oscillations and g-mode ones appears in their prop- 
agation regions. The propagation region of g-mode oscillations with frequency u and azimuthal 
wavenumber m is (u — mQ)^ < k^. That is, the evanescent regions of the inertial-acoustic os- 
cillations (the shaded areas in figures 1, 2 and 3) are now the propagation regions of g-mode 
oscillations. That is, the shaded region in the figure is now the propagation region, and thus 

^ So far, we less discussed the axisymmetric oscillations by assuming that they will not produce much light 
variation. 
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Fig. 1. Frequency-radius relation, uj = (m ± l)f2±K (four combinations of ± being possible), showing 
resonant radii where inertial-acoustic oscillations with frequency u) and azimuthal wavenumber m = 1 
have resonant interaction with disk deformation (thick solid, thick dashed, and less thick dashed curves). 
Horizontal resonance of inertial-acoustic oscillations is considered. In some part of the frequency-radius 
curve of a; = [m ± 1) ± k, the resonance excites inertial-acoustic oscillations, while in the other part the 
resonance dampens the oscillations. The part of excitation is shown by thick solid (or thick dashed) curves, 
while the part of damping is shown by a less thick dashed curve. Even in cases of excitation of oscillations, 
there are cases where the excitation occur in the evanescent region of the oscillations. Such cases are shown 
by thick but dashed curves, since the interaction between the oscillations and the deformation of disks is 
weak there. The evanescent region of m = 1 inertial-acoustic oscillations is shown as shaded area. For 
comparison, the curves of ^{r) and 2f2(r) are also shown. The horizontal dashed curves demonstrate the 
propagation of typical oscillations that are excited. The open circles on the horizontal dashed curves show 
the radii where resonance occurs. The propagation of oscillations labeled by wll is bounded in a finite 
range: The inner boundary is the inner edge of the disk and the outer boundary is the barrier of — k. In 
the case of oscillations labeled by wh, the propagation region is bounded only in the inner boundary. In 
this sense, the latter oscillations labeled by tjfj will be less observable with a definite frequency, compared 
with the former oscillations of t^LL- 

the thick sohd curve and the thick dashed one are interchanged. 

This consideration shows that g-mode oscillations are trapped in a region with a finite 
radial extend, and in their propagation region, a resonance occurs so as to amplify the oscilla- 
tions (see figures 1, 2 and 3). In the case of non-axisymmetric g-mode oscillations, however, we 
should remember the following situations. Except for the cases where the frequencies of g-mode 
oscillations are as high as the upper value limited by the propagation region of the oscillations, 
the corotation radius (c^; = mVt) appears in the propagation region of the oscillations (see fig- 
ures 1 and 2). Different from the case of inertial-acoustic oscillations of n = 0, the presence of 
corotation resonance in the propagation region strongly dampens the g-mode oscillations (Kato 
2003; Li et al. 2003). This strong damping suppresses the growth by the resonance by disk 
deformation considered here. In a narrow frequency range where the frequencies are as high 
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Fig. 2. The same as figure 1 except for m = 2. Because of tlie difference of m, tfie frequency-radius 
relation of resonance is now = f2 ± k and lo = 3f2 ± n (both signs of ± being possible) . On the thick solid 
parts of the curves, resonance excites oscillations in their propagation region, but on the thick dashed part 
resonance occurs in the evanescent region of oscillations. The less thick dashed parts of the curves show 
the resonance that dampens inertial-acoustic oscillations. The evanescent region of the oscillations are 
now specified by 217 — k < a; < 257 + /t (dark area) . For comparison, the curves of fl, 2U, 3fl are also shown. 
The propagation region of oscillations labeled by is bounded in a finite region. However, the region of 
oscillations labeled by t^HH is not bounded in the outer region. In this sense, the latter oscillations will be 
less observable with a definite frequency, compared with the former oscillations with cui^. 

as the upper limit, the g-mode oscillations can be trapped in a radially narrow region without 
having corotation resonance in their propagation region (cf. figure 1). In this case, however, a 
resonance does not occur inside the propagation region (see figure 1). 

To conclude, excitation (amplification) of global non-axisymmetric g-mode oscillations 
is less promising compared with that of global non-axisymmetric inertial-acoustic oscillations. 
Axisymmetric g-mode oscillations are, however, a good candidate of one of observed quasi- 
periodic oscillations. Axisymmetric g-mode oscillations are trapped around ^max and excited, 
since a resonance that can excite the oscillations appears in the propagation region of the 
oscillations (see figure 3). This will be one of possible excitation mechanisms of the stable 67 
Hz oscillation observed in GRS 1915+105. 

8. Discussion 

In this paper we have examined resonant excitation (amplification) of disk oscillations in 
deformed disks. Analyses are made in a general way as far as possible so that we can perspective 
see the essence of the resonant amplification processes. 

In a previous paper (Kato 2004) we have considered a warp as a possible disk deformation 
to excite oscillations. In this paper, we show that, in addition to warps, a one-armed pattern 
symmetric with respect to the equatorial plane can also excite (amplify) disk oscillations. 
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Fig. 3. The same as figure 1 except for m = 0. The dark area represents the evanescent region of 
inertial-acoustic oscillations (the propagation region of gravity oscillations). Oscillations with frequency 
of Lo have resonance with disks at the radius where lo — ^zLkis satisfied. This frequency-radius relation is 
shown by separating into three parts; thick solid, thick dahsed, and less thick dashed parts. On the thick 
solid and thick dashed parts, resonance excites inertial-acoustic oscillations in their propagations region 
and in their evanescent region, respectively. (In the case of g-mode oscillations, the meaning of the thick 
solid part and the thick dashed one is exchanged.) On the less thick dashed part, resonance dampens 
inertial-acoustic oscillations. It is noted that trapped g-mode oscillations whose frequencies are around 
'*max can be excited by our present resonant model. 

As shown in the text, the horizontal resonance is of importance, since it can excite 
inertial-acoustic and g-mode oscillations. The resonance occurs at the radii of the resonant 
condition, [cj — (m± = k^, being satisfied, where uj and m are frequency and azimuthal 
wavenumber of inertial-acoustic (or g-mode) oscillations, respectively. 

Important quantities determining the stability are, (i) the sign of work done by deformed 
disks to disk oscillations through resonance at resonant radius, and (ii) the sign of wave energy. 
As discussed in the text, the sign of the work done on the oscillations with uj and m is determined 
by the sign of — (m -f 1)^2 [oi uj — [m. — 1)VL\ at the resonant radius, depending that the path of 
coupling ism— >m-|-l — 1 ^ m. The sign of wave energy is determined where 

the corotation radius exists. If a wave is mainly inside the corotation radius oi uj — rriQ = 0, 
the wave energy is negative, while it is positive when the main part of the wave is outside the 
corotation radius. By combination of the sign of the work done on waves and the sigh of the 
wave energy, excitation or damping of the waves is determined. 

Whether resonance appears in the propagation region of oscillations or not is also impor- 
tant for observability of the oscillations. The propagation region of inertial-acoustic oscillations 
is specified by {uj — rriQ)'^ > k^, while that of the g-mode oscillations is {uj — mf2)^ < k^. 

In this paper, some of the above results [i.e., the sign of work done on oscillations by 
resonance] are derived under the approximation that the oscillations are local in the sense that 
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their radial wavelength is shorter than the characteristic radial scale of disks (see section 5). 
The results, however, are free from the approximation, and valid even in the case of fcr < 1. 

The results mentioned so far restrict frequency ranges of resonantly excited oscillations, 
but their frequencies are not determined discretely or uniquely. Further restrictions concerning 
excited oscillations are made in section 6 by introducing assumptions that the oscillations 
are local and thus they are predominantly at the boundary of their propagation region, i.e., 
{uj — mVlY ~ /t^. Results based on this assumption are summarized in tables 1 and 2 for two 
cases where the disk deformation is a warp (table 1) and it is a one-armed one plane- symmetric 
with respect to the disk plane (table 2). As shown in these tables, inertial-acoustic and g-mode 
oscillations are excited at the radius of 4rg. Applications of these results to observed high- 
frequency QPOs are made, for example, by Kato and Fukue (2006), Kato (2007) and papers 
referred therein. 

As mentioned before, unless the local approximation of oscillations is introduced, res- 
onantly excited oscillations are not determined uniquely. What is determined is possible fre- 
quency range of excited oscillations. This issue is discussed in section 7 and the results are 
summarized in figures 1, 2 and 3. Eigenfunctions of global trapped oscillations will determine 
the most probable frequencies of resonantly excited oscillations. This problem, however, is 
outside the scope of this paper. 

Finally, it is important to note that for the resonance considered in this paper to occur, 
the disk must be relativistic, i.e., otherwise, the resonant condition is not satisfied in the disk. 
In other words, for resonance to occur in relativistic disks, the disk deformation must be one- 
armed (i.e., a warp or one-armed pattern symmetric with respect to the equator). 

Appendix A General Form of Nonlinear Coupling Terms 

In general cases of Fi ^ 1, the pressure coupling terms, Cp, given by equation (6) should 
be modified, although the term of is unchanged. From equations (9)-(14) of Kato (2004), 
we have, after some manipulations, 

d ( di,di,\ , d 



PoC„.(tO=-^(po^^j + 
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1 d 



2dn 



Fi(Fi - l)po 



(82) 



dr j dvi . 

Then, among and such commnesurable relations as 

(83) 

hold, where ^^^^ and are arbitrary functions of r. This can be shown by performing 
the integrations by parts, so that the integrands become symmetric with respect to ^^^\ ^^'^\ 
and This commensurable relation is a general nature of conservative systems. In the 

cases where varous perturbations are on a disk, energy exchange among these perturbations 
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may occur through resonant processes. The above commensurable relation, however, show that 
even in such cases, the total energy of perturbations is conserved in a conservative system. In 
other words, there is no spontaneous growth of perturbations in equilibrium disks. This is a 
reason why we consider deformed disks as the case of excitation of oscillations. 
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